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ABSTRACT 

In the present study, batch adsorption of dissolved particles in coal effluent was examined. Analysis of 
the process was studied under 3-level full factorial design, with reference to three (3) factors viz: Mass 
of adsorbent, Xj(g), Contact time, x 2 (mins), Volume of effluent, x 3 (ml). In other words, 3 -full factorial 
experiment was carried out. The average responses generated at various random experiments varied 
between 0.0525 - 0.065mg/l, and used to generate a linear model to fit the experimental data. The 
design expert computer software (7.0 version) was used to plot 3-dimensional surface relationship of 
the factors. It was observed, from the interaction plot of effluent volume versus adsorbent mass, that as 
the mass of adsorbent increases (with decrease in the volume of effluent), the uptake of the total solid 
particles rises. The case is different for the interaction of adsorbent mass and contact time, where more 
particles are being removed as both mass of adsorbent and contact time are increased. The variance 
(S u 2 ) for the various number runs fall between - 4.5xl0" 6 . The Chochrain's distribution (G) was 
calculated to be 0.3871. By simple comparison of G ca kuiated with G ta bie (0.6798) at 0.05 level of 
significance, it was deduced that the hypothesis for the homogeneity of variance is acceptable. The t- 
test was, also, used to confirm the viability of the coefficients of the model equation, while the F-test 
was used to assess the adequacy of the model. The value of F ca i cu i ated was found to be 0.7646, while that 
of F ta bie was 4.07 (>F ca icuiated), and this confirmed the adequacy of the model. 
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INTRODUCTION 

Design of Experiment (DOE) is the design of all information-gathering exercises, where variation is present, 
whether under the full or partial (fractional) control of the experiment. In testing, each experiment gets a 
minimum amount of measurable conversions, known as the sample size per data, for any statistical analysis 
(Billing, 2008). Therefore, the more experiment you have, the more conversions you require. The conversion 
data can also be thought as time, since the longer you leave your web page up, the more data you get. In an 
experiment, we deliberately change one or more process variables (or factors) in order to observe the effect 
(change) on one or more response variables (dependent variables). The design of experiments is an efficient 
procedure for planning experiments so that the data obtained can be analyzed to yield valid and objective 
conclusions. Hunter and Hunter (1998) reported that DOE begins with determining the Objectives of an 
experiment and selecting the process factors for the study. According to the report, an Experimental Design is 
the laying out of a detailed experimental plan in advance of performing the experiment. Well chosen 
experimental designs minimize the amount of information that can be obtained for a given amount of 
experimental effort. The Factorial Experiment is, thus, an experiment whose design consists of two or more 
factors, each with discrete possible values or levels, and whose experimental units take-on all possible 
combinations of the levels across all the factors (Maiti, 2001; Kashik, 2007). Brathy (1998) and Shilpi (2003), 
also, reported that a factorial experiment allows the study of the effect of each factor on the response variable, as 
well as the effects of the interactions between the factors on the response variable. During any experimental 
design, all possible interactions may be included in the test (Full Factorial), or some of the interactions are 
selectively used (Fractional Factorial). 

For the vast majority of factorial experiments (Dick and Ewing, 1997; John, 2001; Box and Wilson, 1951), each 
factor is assigned two levels-f/ie high and low levels. But for computational purposes, the factors may be scaled 
such that the high level is assigned a value of +1, and the low level is assigned a value of -1. In line with this 
arrangement, a full factorial design contains all possible combinations. Since full factorial gathers additional 
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information, it reveals all possible interactions, but there is always a "trade-off when a higher number of factors 
is involved (>5 factors); more data equals more information, but implies longer test duration with very high data 
requirements, since it shows all the experiments. On the other hand, the fractional factorial includes everything 
from only main effects to the interactions. 

Frank Yates made significant contributions, particularly in the analysis of designs, by the "Yates Analysis" 
(Wikipedia, 1990). The Yates analysis exploits the special structure of the designs to generate least square 
estimates for factor effects, for all factors and relevant interactions. Maiti and Kenway (2004) reviewed that the 
analysis can be used to determine the ranked list of the factors, as well as the goodness-of-fit (as measured by 
the residual standard deviation) for the various models. The document unveiled that any given data should be 
arranged in "Yates order" before performing the Yates analysis. In this regard, given 'K-factors', the 'Kth' 
column consists of 2 (kl) minus sign (that is the low level of the factor), followed by 2 (kl) plus sign (that is the 
high level of the factor). For instance, for a full factorial design with three factors, the design matrix is as 
follows: 



+ - - 

+ 

+ + - 
+ 

+ - + 

+ + 

+ + + 

And the Yates analysis generates the following output: 

A factor identifier (from Yates order)- the specific identifier will vary, depending on the program used to 
generate the Yates analysis. Data plot, for instance, uses the following for a 3-factor model: 

1 for factor 1 

2 for factor 2 

3 for factor 3 

12 for interaction of factors 1 and 2 

13 for interaction of factors 1 and 3 
23 for interaction of factors 2 and 3 
123 for interaction of factors 1, 2 and 3. 

A rank list of important factors - this is the least square estimated factor effects ordered from largest in 
magnitude (most significant) to the smallest in magnitude (least significant). 

A t-value for the individual factor effect estimates, which is computed as follows: T = e/S e 

Where e = the estimated factor effect 

S e = the standard deviation of the estimated factor effect. 

The residual standard deviation that results from the model with the single term only, which is given by: 

Response = Constant X 0.5x,Where X; = the estimate of the /th factor or interaction effect. 

The cumulative residual standard deviation that results from the model using the current terms plus all the terms 
preceding the current term, and is given by: 

Response = Constant + 0.5(all effect estimate down to t, including t). 

This consists of a monotonically decreasing set of residual standard deviations (indicating a better fit as the 
number of terms in the model increases). The first cumulative residual standard deviation (response = constant) 
is for the model, and the constant is the overall mean of the response variable. The last cumulative residual 
standard deviation is given by: 
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Response = constant + 0.5 (all factor and interaction estimates). 

It is important to note that if there are k-factors, each at 2 levels (high and low), then a full factorial design has 
2 k runs (Perry et al, 1984; Sinnott, 1999). Following from this, when the number of factor is 5 or greater, a full 
factorial design requires a large number of runs, and is not very efficient. In this direction, the design guideline 
table recommends a functional factorial design or a plackett-Burman design for such case (Rose and Garba, 
1991). The research work in course, thus, explores the content of this recommendation in the assessment of 
adsorption capacity of feather-derived adsorbent in wastewater management applications. 

MATERIALS AND METHODS 
Material Collection/Preparation 

The wastewater sample (coal effluent) was collected from the effluent holding pond of the Coal Mining 
Industry, Akwukwe-Enugu, located at the South-Eastern part of Nigeria. The soil in the region is rocky with 
high coal deposit. Concentration of particulate matters in the coal washery effluent is always high as a result of 
its exposition to heavy down pours of rain. 

The feathers of a Chicken were picked from the chicken-slaughter points at Obinze market in Imo State of 
Nigeria, sliced and impregnated in 60% strength of orthophospheric acid for 24-hours for activation. At the 
elapse of the time, the feather sample was removed and washed severally with distilled water to obtain a sample 
of neutral or nearly neutral pH (using a Delta-320 pH meter). The sample was then inoculated in clean trays and 
dried by means of oven dryer. It was, further, inoculated in clay pots and subjected to carbonization process at a 
charring temperature of 300 C C, using a furnace equipment of model- 2KX0-60. A ceramic mortar set was used to 
grind the adsorbent sample (which is now in the form activated carbon), and a sieve size of 600{im was used to 
obtain the granulated size that was used in the study. 

The Factorial Experiment 

Batch adsorption experiments were conducted randomly at three different levels - the High level (+1), Low level 
(-1) and Base level (0), with reference to three experimental factors - Mass of adsorbent (x t ), Shaking/Contact 
time (x 2 ) and Volume of effluent (x 3 ). Initial effluent concentration was checked by means of Ultraviolet 
Spectrophotometer (UVS) before the process. The finally concentration was, also, examined after every random 
treatment, and this served as response of the design process. Eight different experiments were conducted using 
the design array shown in Table 2 (Wikipedia, 1990), with replications to ascertain homogeneity of the 
experimental method. 

RESULTS AND DISCUSION 

The design expert computer software (7.0 version) was used to plot the 3-dimensional surface relationship of the 
factors - Effluent volume (ml), Adsorbent mass (g) and Contact time (mins.). It could be observed from the 
interaction plot of effluent volume versus adsorbent mass (Fig 1) that as the mass of adsorbent increases (with 
decrease in the volume of effluent), there is higher uptake of the total solid particles, unlike in figure 2, where 
more particles are being removed as both contact time and adsorbent mass increase. Figure 3 follow the same 
trend as Figure 1, that is as time increases (with decrease in effluent volume), there is more uptake of the 
particles. This suggests that low effluent volume implies less particle concentration in the solution. 

Table 1: BASIS FOR THE FACTORIAL DESIGN 



Factors 


High Level (+) 


Low Level (-1) 


Base Level (0) 


Mass of Adsorbent, x^g) 
Shaken time, X 2 (min) 
Volume of Effluent, X 3 (ml) 


1.0 

10.0 
20.0 


0.2 
2.0 
10.0 


0.6 

5.0 
15.0 



The proposed model equation is given by: 

Y u =bo+b 1 X 1 b2X2+b3X 3 +bi2XiX 2 b 1 3X 1 X3+b23X2X3 



For each of the experiments (with replication, i.e. n=2), check for mean value of the runs (y), the variance (S n ) 
and the Chochrain's distribution (G); 
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Where: y u = (yi+yi)/2 

S u 2 = ' Z (v„, -v„) 2 and 
G = Max. S - 2 ! jj? =1 

u=l 
These may be arranged in a summary table of array as follows: 

Table 2: Design Array for Factor Responses 



N 


fn 


f, f, 


f, £, U 


ffi 






Y^mg/l) 


Y 2 (mg/1) 


^l(mg/l) 


s u 2 


An 


Xl X2 


X.3 X.!X 2 XjXj 


x 2 x 3 


1. 


+1 


+ 1 


+ 1 +1 


+1 


+1 


+1 


0.055 


0.056 


0.0555 


2xl0" 6 


2. 


+1 


-1 


+ 1 +1 


-1 


-1 


+1 


0.050 


0.055 


0.0525 





3. 


+1 


+ 1 


-1 +1 


-1 


+1 


-1 


0.056 


0.055 


0.0555 


4.5xl0" 6 


4. 


+1 


-1 


-1 +1 


+1 


-1 


-1 


0.050 


0.056 


0.064 


2x1 0" s 


5. 


+1 


+ 1 


+ 1 -1 


+1 


-1 


-1 


0.055 


0.054 


0.065 





6. 


+1 


-1 


+ 1 -1 


-1 


+1 


-1 


0.051 


0.055 


0.064 


5xl0 -7 


7. 


+1 


+ 1 


-1 -1 


-1 


-1 


+1 


0.055 


0.053 


0.065 





8. 


+1 


-1 


-1 -1 


+1 


+1 


+1 


0.052 


0.055 

VTr. 1 A 


0.0635 


4.5xl0" 6 



Thus, G : 



4.5x10'' 



: 0.3871 

1.1625xl0" : 



But G table (0.05,8,1) = 0.6798 

Therefore, since G Tab | e = 0.6798>G ca | CU i ated = 0.3871, the hypothesis for the homogeneity of variance is 
acceptable. 

NB : The variance for the individual experiments is given by: 
S E 2 = _1_ ZS U _1_ (4.65X10 -5 ) = 5.8125X10 6 



Evaluating the coefficients of the model equation: 
b„= 1 (0.4315) = 0.0539 

8 
b, = _1_ (0.0555-0.0525+0.0555-0.053+0.0545-0.053+0.054- 

8 0.0535) =9.375xl0 3 

b 2 = _!_ (0.0555-0.0525+0.0555-0.053+0.0545-0.053+0.054- 
8 0.0535) =6.25xl0- 5 

b 3 = _1_ (0.0555+0.0525+0.0555+0.053-0.0545-0.053-0.054- 
8 0.0535) =6.25xl0- 5 

bi2 = _1_ (0.0555-0.0525-0.0555-0.053+0.0545+0.053-0.053+ 
8 0.0535) =1.875xl0' 3 

bi 3 = _1_ (0.0555-0.0525+0.0555-0.053-0.0545+0.053-0.054+ 
8 0.0535) =4.375xl0' 3 



J_ (0.0555+0.0525-0.0555-0.053-0.0545-0.053-0.054+ 
8 0.0535) =-6.255x10'' 

t tab ie = (8,0.05) = 2.3061 
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■N 



5.8125x10'' 



■■ (3.6328) 1 ' 2 = 6.0273xl0' 4 



16 



Thus, (table) . (S b ) 



: (2.3061) (6.0273xl0" 4 ) = 139x10 s 



Comparing the value with coefficients of the model equation: 

(ttabie) ■ (Sb) = 1.39xl0" 3 >b 2 , bi2 and b 23 ; hence, they dropped from the model. 

The model equation then becomes: 

Y„ = bo+bjXj + b 3 X3+b 12 x 1 X2+bi3X 1 X3 

Y u =0.0539+0.009375x 1 +0.001875x 3 +0.991875x,x 2 +0.004375x 1 x 3 

Table 3: Checking for the Adequacy of the Model 



N 


l 


2 


3 


4 


5 


6 


7 


8 


y u 


0.0714 


0.0408 


0.0677 


0.0439 


0.0598 


0.0545 


0.0552 


0.0489 


y; 


0.0555 


0.0525 


0.0555 


0.0530 


0.0545 


0.0530 


0.0540 


0.0535 


y u -yu 


0.0159 


-0.0117 


0.0122 


-0.00091 


0.0044 


0.0015 


-0.0012 


-0.0046 


(y» yu) 2 


2.53xl0" 6 


2.53xl0" 4 


1.37xl0" 4 


2.28xl0" 5 


1.94x10 s 


2.25xl0" 5 


1.44xl0" 6 


2.116x10 s 



SAMPLE CALCULATIONS 

y u ,=0.0539+0.009375+0.001875+0.001875+0.04375 = 0.0774 
y u2 = 0.0539+0.009375+0.001875-0.001875-0.04375 = 0.0408 
y u3 = 0.0539+0.009375+0.001875+0.001875+0.004375 = 0.0677 
y u4 = 0.0539-0.009375+0.001875+0.001875-0.04375 = 0.0439 
y u5 = 0.0539+0.009375-0.001875+0.001875-0.04375 = 0.0589 
y u6 = 0.0539-0.009375-0.001875-0.001875+0.004375 = 0.0545 
y u7 = 0.0539+0.009375-0.001875-0.001875-0.04375 = 0.0552 
y u8 = 0.0539-0.009375-0.001875+0.001875+0.04375 = 0.0489 
Then, I(y„ - y u ) 2 = 666605x10'" (from the table). 
Thus, Q L =2x6.66605xl0" 6 = 1. 3321x10'' 
And S L 2 = Q, 



-1.3321x10- 



Tl 



Thus, 



4,444x10'' 



4.444x10" 
3 

= 0.7646 



5.8125x10"' 



But F table (0.05, 8, 3) = 4.07 

Since F=0.7646<F tab i e =4.07, the model is said to be adequate. 
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Volume of Effluent(ml) 



Mass of Adsorbent(g) 



Fig 1: 3-D Surface Plot for Effluent Volume versus Adsorbent Mass 




D.S . 05 

Time(min) 



Mass of Adsorbent(g) 



Fig 2: 3-D Surface Plot for Adsorbent Mass versus Time 




Volume of Effluent (ml) 



Time(min) 



Fig 3: 3-D Surface Plot for Effluent Volume versus Time 



CONCLUSION 

Once the primary variables (factors) that affect the responses of interest are known, a number of additional 

objectives may be pursued. These include maximizing or minimizing a response, reducing variation and seeking 
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multiple goals. What each of these purposes has in common is that experimentation is used to fit a model that 
may permit a rough approximation to the actual requirement. In the present study, for instance, optimum 
average response (0.065mg/l) was recorded at runs 5 and 7; and the empirical model fitted the experimental data 
linearly. 
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